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Abstract 

The notion of Laplace invariants is transferred to the lattices and dis- 
crete equations which are difference analogs of hyperbolic PDE's with 
.«— ■ s ■ two independent variables. The sequence of Laplace invariants satisfy the 

^vj ' discrete analog of twodimensional Toda lattice. The terminating of this 

f**^ ■ sequence by zeroes is proved to be the necessary condition for existence 

Q\ ' of the integrals of the equation under consideration. The formulae are 

^\ _ presented for the higher symmetries of the equations possessing integrals. 

The general theory is illustrated by examples of difference analogs of Li- 
ouville equation. 



1 Introduction 



Liouville equation and its generalizations are studied about one hundred fifty 
years and very many results, including classification ones, are obtained in this 
^ ■ field [jy - |L4| . Surprisingly, it seems that in the discrete version these investiga- 

J-j \ tions were not undertaken, although many results can be transferred from the 

continuous case without any efforts. The aim of our paper is to advance in this 
direction. The equations with explicit dependence on the discrete variable are 
most difficult and for the present we consider only shift-invariant equations al- 
though this condition is, probably, too restrictive. Overcoming of this difficulty 
as well as deduction of classification theorems require further researches. 

Let us consider differential-difference and totally discrete equations of the 
form 

Ui+l, x = f(x, Ui + l,Ui, Uj jX ), (1) 

Ui+1,3 + 1 =/(Ut+l,j,t*i,j+l,'Ui ) j)- ( 2 ) 

Further on we assume that for the equation (pi) f Ui ^ 0, and the right hand 
side of the equation (0) depends explicitly on each of the indicated variables. 



Apparently these equations can be treated as difference approximations of the 
equations of the form 

u xy = f(x,y,u,u x ,u v ). (3) 

They arise also as Backlund transformations and nonlinear superposition prin- 
ciples, and of course this concern both generalizations of the Liouville equation 
and equations integrable by Inverse Scattering Transform (the well known ex- 
amples are dressing chain pi and difference KdV equation [p7[). 

One can divide the different generalizations of the Liouville equation into 
several (intersecting) subclasses, subject to what property of the Liouville equa- 
tion is chosen as definition of integrability, for example: explicit formula for 
the general solution, differential substitution to linear equation, symmetries de- 
pending on arbitrary functions, or integrals along both characteristics. In the 
present paper we accept the later property as definition. As far as we know 
it was offered as independent definition in the paper M and the corresponding 
subclass was called 'Liouville type equations'. Since this term is suitable rather 
to join (or intersection) of all subclasses mentioned above, we would say just 
'equations possessing integrals'. Other term, 'Darboux integrable equations' 
[ fill , [14 [l3| is not exact, since Darboux had used rather the first of the listed 
properties. 

The given definition means that the lattice (Q) admits the x-integral X 
and i-integral /, that is functions on x and finite set of dynamical variables 
{u i+m ,u\ n ' — d"(ui)\ m G Z, n e Z + } (i is assumed to be fixed here), such 
that the identities hold in virtue of the lattice 

D x (X) = 0, (Ti-l)(I) = 0. 

The sequence of the coupled Riccati equations 

Ut+i, x - Ui, x = uf +l - u\ (4) 

can be considered as the classical example, its integrals are of the form 

T 2 Y 1 n (uj -u l+2 ){u l+ i -u l+3 ) 

1 = Uix - U i: X = {Ui,Ui + i,Ui + 2,Ui + 3) = ~, r~. r. 

[Ui ~ u i+3 )(u l+1 - u l+2 ) 
Slightly more complicated example is given by the lattice 

Ui+l, x Ui, x = Ui+i + u t (5) 

for which, as one can easily prove, 

2/2 v _ ( M *+3 - U i+ l)(u i+ 2 ~ Ui) 



1 = (u i>xx - i)7<», x 



Ui+2 + Ui+1 



It is not difficult to prove that in general case / depends only on derivatives 
u\ n ' , while X depends only on shifts u i+m . We will call the order of the higher 



derivative on which the i-integral depends the order of i-integral, and the order 
of x-integral is equal to the difference between the maximal and minimal shifts 
of the variable m which are involved in it. So, in the last example ordl = 2, 
ordX = 3. 

To solve the lattice possessing integrals one has to solve the difference equa- 
tion X = Ci where Cj is an arbitrary sequence, and consistent with it ordinary 
differential equation / = C{x) where C(x) is an arbitrary function. Usually this 
allows to obtain the explicit solution as the final result, for example, one has for 
the lattice (§) 

YXX , tyx 

2ip x Zi - tp 
where function ip{x) and parameters Zi are arbitrary. Notice, that 

nr _ n _ YXXX &W xx V / \ 

11 — b^ — — - 9"72~' ~~ ^ Zi ' Zi +!' Z i+2; z i+3>- 

Yx ^Vx 

Nevertheless, the relation between existence of integrals and explicit formula for 
solution remains unclear till now even in the most studied continuous case. 

For the discrete equation (|2]) our definition of integrability means existence 
of i- and j-integrals satisfying relations 

(T, - 1)(7) = 0, (T, - 1)(J) = 0. 

Now the dynamical variables are Ui +m j, Uij +n , m,n € Z, since the rest can 
be eliminated in virtue of (0). ft is easy to prove that z-integral depends on 
the shifts Uij+ n only, and j-integral depends on the shifts Ui+ m ,j. The order of 
integrals is defined in just the same way as order of x-integral for the equation 
(Qj). As an example, let us consider the equation 

Ui+i,j+i - Uij + i -\ = (6) 

which admits z-integral of the first order and j-integral of the third order: 

j_ ! T _ (Ui+M ~ u i+2,j)(ui+l,j - Uij) 

1 — u i,j+l H 1 J — 7 T7 7- 

Uij {Ui+3,j - Ui + l t j)(Ui + 2,j ~ Uij) 

Now we summarize our results. In the Section 2 the notion of the Laplace 
invariants is transferred to the discrete case (cf. iQ). The sequences of Laplace 
invariants for the above examples and discrete Liouville equations from the Sec- 
tion 3 are terminated by zero and this allows us to conjecture, by analogy with 
continuous case, that this property is characteristic for the equations possessing 
integrals IO] -- [H. The main result of our paper is the proof of this conjecture 
in one direction (necessity), given in the Section 5. The outline of the proof 
is the same as in continuous case. In Section 6 the constructive formulae are 
presented for the symmetries of the lattices (H) possessing integrals. 



Laplace invariants satisfy some recurrent relations analogous to the discrete 
Toda lattice p8| - p3fl . In Section 4 these relations are identified as the well 
known formulae for the Backlund transformations and nonlinear superposition 
principle for the two-dimensional Toda lattice. The proof of sufficiency in the 
continuous case is based essentially on the results on the structure of the inte- 
grals of this lattice ra, [L5| and therefore it would be very important to obtain 
analogous results in the discrete case as well. 

2 Discrete Laplace invariants 

Laplace invariants of the linearized equation are very useful tool in the theory 
of the equations (||). In particular, it was proved that the equation possesses 
integrals iff the sequence of the Laplace invariants hj is terminated by zero at 
the both ends Q Q. The invariants h^i,h are calculated directly by the 
right hand side of the equation and the rest are found recurrently by the formula 

(\oghj) xy = hj+i - 2hj - + hj-i (7) 

which is equivalent to the two-dimensional Toda lattice. 

It turns out that analogous constructions can be defined in the discrete 
situation as well. Firstly let us consider differential-difference operator 

L = T i D x + a i D x + b i T i + c i = (T< + ai )(D x + ^_i) + H ifi = 

= {D x + bi)(T i + a i )+Ki,o 

where 

Hi.o = Ci - ai&i_i, K h0 = a- a l b l - D x (ai). 

One can easily check that if Hi$ ^ then the formula holds 

{D x + b i , 1 )L = L 1 (D x + b i _ 1 ) (8) 

where 

Li = (D x + bi,i){Ti + a,) + H i>0 , bi A = 6 4 _i - D x (logFi, ). 

It means that the operator D x + bi^\ maps the kernel of L into the kernel of the 
operator L\. Moreover, operator -jj—iTi + a^) define the inverse transformation 
which maps the kernel of L\ into the kernel of L : 

L——(Ti + a,i) = (Ti + a l )——L 1 . 

We shall call the mapping L = L ^ Li the Laplace x-transformation. Its 
iterations yield the sequence of operators 

L 3 = (T t + ai )(D x + bi-u) + H ij: j > (9) 



where Hij, bij are defined by recurrent formulae 

Hi,i+i — Hij + D x {a,i) + ai{bij + i — 6j_ lj+1 ), 

b iij+ i = k-ij - D x (\ogHij), b it o = h. (10) 

Laplace ^-transformation is defined in the similar manner. Analog of the 
formula (||) looks 

(Ti + cn,i)L = L-i(Ti + at) (11) 

where 

L-i = (Ti + a. h i)(D x + bi) + -FQ+1,0, 0»,i = a^i+1,0/^,0 

and it is assumed that liT^o ¥" 0- Operator -^—(D x + bi) defines the inverse 
transformation: 

L-—(D x + bi) = (D x ' 



Kifi Ki+1,0 

Iteration of the ^-transformation generates the sequence of the operators 

L-i = (D x + b i+j ) (Tt + a l3 ) + Kij , j > (12) 

where K^, a%j are defined by recurrent relations 

Ki,j+i = &i+ij ~ D x (dij+i) — di,j+i (6i+j+i — bi + j), 

o-Lj+i = OijKi+ij/Kij, a-ifi — ai. (13) 

It is more convenient for calculations to consider, under assumption a% 7^ 
(this guarantees o»j 7^ 0), the values 

h.. — I Hij/ai, j > 0, 

-j-l/ a i,-j 



\ Ki_j_x/ai-j-i, j<0 



which we will call Laplace invariants. Actually, this exchange is necessary, 
since the operators L should be considered up to the gauge L = a~ +1 Lai and 
the values hij are just the invariants of this transformation. In particular, 
comparison of the formulae (|3|) , (|ll| ) and their inverses demonstrates that the 
Laplace x- and ^-transformations are inverse to each other up to the gauge 
transformation. One can easily check that hij satisfy the lattice 

log n ) = h,j-i - h^ - hi+ij + h i+1J+1 (14) 

Ki + i,j J x 

while the values Hij, K^ do not satisfy closed system of equations. 

Laplace invariants of the lattice (0) are understood as Laplace invariants of 
the linearization operator 

L* -M-^Ac JUi x X /(Ijll^i JUi ■ 



Passage from the values Hij and Kij to hij is always possible in virtue of the 
nondegeneracy condition f Ui x ^ 0. Notice also, that the property ftj+ij = 
Ti(hij) holds due to the shift-invariance of the lattice and therefore the Laplace 
invariants can become zero only on the whole line j = jo- 
Let us return to the examples presented in Introduction. One can easily 
check by direct calculation that for the lattice (|4|) the invariants hi_ 3 and /i^o 
become zero, as well as invariants hi_ 3 and hii for the lattice (g). So the 
conjecture arises that the criterion of the existence of the integrals mentioned 
above is valid also in the discrete case. 

Laplace transformations and invariants for the difference-difference operators 
arc defined quite similarly. Moreover they were studied in the review |y| and 
we only bring few key formulae. Let us consider the operator 

L = TiTj + ciijTj + bijTi + cij = (Tj + a,ij)(Tj + bi-xj) + -ffjj,o — 

= (Tj + bij)(Ti + Oij-i) + Kij t o 

where 

-^zj',0 = C-ij Of%j0i—X,j^ J^ij,0 = Cij Q>i,j — lVij ■ 

Operator Tj + &j-i,j maps the kernel of L into the kernel of the operator 

L\ = {Tj + b,iji){Ti + dij) + i?ij+i,0) hj,i = bi-ijHij+io/HijQ 

and operator j, 1 (T + a^) defines inverse transformation. Iterations of this 
transformation, which we will call Laplace j-transformation, yield the sequence 
of operators 

Lk = (T + ai,j+k)(Tj + h-ijk) + Hijk, k > 
where i?^, biju are defined recurrently: 

Hij.k+l — Hi,j + l,k + &i,j+kbij,k+l — O'iJ+k+lbi-lJ.k+l, 

bij : k+i = bi-\_jk.Tli.j+\,ki Hijk, bij t Q = bij. 

Laplace i-transformation and values Kijk, aijk are defined analogously. 

When dijbij ^ one can introduce Laplace invariants hijk — H^k / ' (ai,j-\-kbi—i,jk) 
for k > 0, hijk = Kij,—k—i/(ai,j—i,—k—ibi—k—i,j) for fc < 0. As in the previous 
case these quantities do not change under the gauge transformations and satisfy 
closed system of equation 

(^j+i,fc-i + l)(^+ij,fc+i + l) = (^ +lfe + l)(V + 1 lj - fe + l)^fe^+ij+i,fc- (15) 

It is easy to check that for the example (|6|) invariants fry,-3 and hij t o of the 
linearized equation become zero in accordance with our conjecture. 



3 Discrete Liouville equations 

Remind that the Liouville equation [pi can be obtained from the Toda lattice 
(0) ^ imposing the simplest boundary conditions h\ = h-\ = 0. The natural 
question arise, what is the result of the same reduction of the discrete Toda 
lattices (|14| ) and ( |l5| ) . We will see that structure and properties of these three 
equations are very similar: each of them possesses second order integrals, has 
only two nonzero Laplace invariants, and admits linearizing substitutions of the 
first order. Moreover, the passage to the limit turns the discrete versions into 
the Liouville equation itself. Let us consider these examples in more details. 

Boundary conditions hi t \ = h^-i = turn the equation ([14]) into the lattice 
of the form (yj) on the variables Ui = hi t o : 

Ui + i tX Ui - u i+ iu iiX = u i+ iUi(u i+ i + m). (16) 

It is easy to check that assuming Ui(x) = eu(x,y), y = ei and passing to the 
limit £ — » one obtains the Liouville equation (logu) xy = 2u. Notice also that 
the lattice (16) can be obtained from (4J) by means of substitution Ui = Uj+i— Uj. 



The linearization operator for (16) is 



L = D x T t — D x - (u lx /v,i + 2u i+ i + Ui)Ti + u i+ i(u itX /uf - 1) 

Ui 

and calculating of the Laplace invariants gives hii = h%.-2 = 0. 
Separating the variables Ui and Ui+\ and differentiating yields 

2 2 

Ui-\-\x V*i,x . V>i-\-l,xx i-\-l,x ^i,xx i,x . 

Ui+1 = + Ui, ^ U l+ i. x = : 2" + U ix 

U l+ i Ui U l+ i Uf +1 Ui uf 

from which the i-integral can be easily found: 

T ~Ui xx ^i,x 2 

1 = 2 3 — ^ m, . 

u t uf 

It is not difficult to find the ir-integral as well: 

X = (1 + Ui/u. l+1 )(l + m/ui-i). 

Equation ( |16| ) is obtained from the linear equation Vi+i yX — Vi,x = by means 

of the substitution 

(v i+ i - Vi)v ix 

Ui = — 

Vi+lVi 

(cf. substitution u ~ v x v y /v 2 which linearize the Liouville equation) and this 
allows at once to write down the formula for the general solution: 

, s {C i+ l - Ci)lp x 

Ui{X) = 



(Cj+l +1p)(Ci + 1p) 



where a are arbitrary constants and ip( x ) is arbitrary function. 

In the totally discrete case the boundary conditions hij t i — hij t -i = reduce 
the equation (15) to equation of the form (0) on the variables Uij — hij t o : 



u i+1 j +1 (l + l/u l+1J )(l + l/ui ij+ i)uij = 1. (17) 

It turns into the lattice (jig ) after the passage to limit Uij — —sui(ej), e — > 0. 

The Laplace invariants of this equation hij t i and /iy,-2 are equal to zero just 
as in the previous example. The integrals are given by formulae 



1 = 1 



Uij(u i} j-l + 1) 



a 



1 j 



Ui,j-l J V u ij+l( u ij + 1) 



J=\ 1 

The substitution 



Uij(Ui-lj + 1) 



u 



v 



Ui-l t j J \ U i+ ij(Uij + 1) 

(Vi+l,j -Vij)(Vi !J+ i -Vij) 



brings ( |l7[) to the linear equation Vi+i j+i —Vi+ij —Vij+i +Vij = and therefore 
the formula for the general solution is valid 

(cj+i - Cj)(k j+1 - kj) 

lJ (C»+1 + fcj)(Ci + fcj+l)' 

where q, fcj are arbitrary constants. 

4 Discrete Toda lattices 

Now we slightly digress from our main subject in order to discuss the equations 
( pT[ ) and ([l5]) which are the natural discrete analogues of the Toda lattice (^) . 
These equations have been already studied in literature from this point of view 
(possibly in somewhat different notation) see e.g. |L8| - [g3|. We should like to 
reproduce for completeness some well-known formulae which demonstrate that 
the links between these three lattices are more close than it seems at the first 
glance. The results of this Section will not be necessary in what follows. 
It is more convenient to start from the lattice 



' 3,xy 



exp(r i+ i-2r i +r i _i), (18) 



rather than from (M) which is related with it by means of substitution hj = Tj xy . 
The Backlund transformation for this lattice is of the form 

(ry - r i+ ij) x = exp(r ij _i - r i3 - r i+ i tj + r l+1J+1 ), (19) 

(n+i,j - n,j-i)y = exp(fij - r i+ i t j - r irj -i + r 4+ i,j-i)- (20) 



Indeed, the cross-differentiation yields the formula 

rij.xy = exp(r ij+ i - 2nj + nj-i) + c(x,y) 

where function c(x, y) does not depend on i,j and can be easily removed by 
shift of the form r — > r + C(x, y). The lattice (119) is reduced to (n4) by means 
of substitution hij = (?*y — ri+i t j) x . 



We see that the second subscript in the formulae Q19J), (20) corresponds 
to the shift in the lattice (O) while the first one denotes the iteration of the 
Backlund transformation. Now introduce the third subscript and consider the 
Backlund transformation which changes it and leaves the first subscript fixed 
(vice versa, the fixed third subscript must be assigned to the variable r in the 
formulae ©, pf 



(rijk ~ rij,k+i)x = exp(ri,j_i,jfe - njk - r^.k+i + nj+i,fc+i), 
(r%j,k+i - n,j-i,k) y = exp(ryfc - n^u+x - r i,j-i,k + n,j-i,k+i)- 

Let us consider the superposition of these Backlund transformations. It is suf- 
ficient to consider only x-part, since differentiating on y brings to the same 
result. Calculating (r^k — rj+ij ] fe+i) x in two different ways one obtains after 
some transformations the relation 

( T j - 1 )( ex P( r i+ij-i,fc - rj+ij.fe+i - ri,j_i,fc + rij,k+i) - 

- exp(n jJ -_i ] fc + i - n+i^k+i - n,j-i,k + n+i,jk)) = o. 

The expression in the brackets is equal to some constant dk which can be set to 
1 by shift of the form njk —* fijk + Cjfcj resulting in desired difference equation 

exp(r,+i,j_i,fc + r ijt k+i) - exp(r,,j_i,fe + i + r i+1 ^ k ) = 

= exp(r, + i,j,fc + i + r»j_i,fc)- (21) 

It is equivalent, up to the linear transformation of the coordinate axes, to the 
well-known Hirota-Miwa equation fill [20l 

Clfi+l,jkfi-l,jk + bfij+l^kfij-l.k + cfij,k+lfij,k-l = 

(it is clear that all coefficients can be removed by means of scaling / — ► a*^* - ^jS 3 '^ -1 )/). 
Substitution 

hijk = exp(r i+ i iiifc +i - ^j-i,fc+i - n+i,jk + n,j-i,k) 
maps equation (K^) into equation 

(h+i,j,k+i + l)(^ij-i.fc + 1) = (K,j-i,k+i + l )(K+i,jk + l)hi+i,j-i,khij,k+i 
which coincides with (HS) up to the linear change of i,j, k. 



It should be stressed that the changes from the variable r to the variable h 
are different in all three cases, and therefore one cannot say that the lattice (ft4j) 
defines the Backlund transformation for (R) and the equation ( |l5| ) is its nonlinear 
superposition principle. Nevertheless, we have demonstrated that these three 
equations live on the same three-dimensional lattice generated by the Backlund 
transformations of the lattice (Jlq). 

5 Necessary conditions for existing of integrals 

In this Section we prove the main theorem which allows to check, rather effec- 
tively, if the given equation (|l|) or (g) possesses the integrals. 

We consider in details only the lattices ([j]). Let g be function on x and 
dynamical variables Ui+ m , u\ n . Denote g* its linearization, that is operator 

oo oo 

.9* = E 9u i+m Tr + Y,9u^ D *- 

m— — OG n—1 

Intermediate calculation proves that the properties hold 

(D x (g))* = D x g* mod [L], (Ti(g))* = T ig * mod [L] (22) 

where L is linearization operator of the lattice (||) and [L] denotes set of all 
operators of the form ML with arbitrary differential-difference operators M. 

Theorem 1 // the lattice (Jip possesses the x-integral X of the order n then 
its Laplace invariant hi —j vanishes for some j, < j < n. Analogously, if it 
admits the i-integral I of the order n then the invariant hij vanishes for some 
j, < j < n. 

Proof. Let us prove by contradiction the first statement. Assume that hi —j ^ 
for all j = 1, . . . ,n. It is the same as the quantities Ky ^ for all j — 
0, . . . ,n — 1 and therefore the formulae (|l2|), (13) define the operators L_j for 



j = 0, . . . , n. The relation (lllf) for these operators takes form 

(Tj + a,i t j)L-j + x = L-j(Ti + a^j-i), j = 1, . ..,n. 

The repeated application of this formula brings to L_jAj_i = mod [L] where 
the difference operators Aj are defined as follows 

A-i = 1, Aj = (Ti + ai^Aj^, j = 0, . . . , n. 

Therefore one obtains the relation 

D x Aj = -b 1+0 A - KtjA^i mod [L], j = 1, . . . ,n. 



10 



Due to the shift T one can assume without loss of generality that X depends 
on Ui, . . . ,Ui+n. Consider the expansion of X* over the operators Aj : X* — 
Y^i=o £?' A7-1 an d multiply it by D x from the left. Collecting the coefficients 
on Aj and D x in the relation D X X* — mod [L] brings to the sequence of 
equations 

(Ac-fci+n-lXCn) = 0, 

(D x - b i+ j-i)(€j) = tj+iKij, j = 1, . . . ,n - 1, 

£o = 

from which the contradiction follows: £o = • • • = 6i = =>• X* = 0. 

The proof of the second statement is quite analogous and we just outline it. 
Assume that ffy ^ for all < j < n and define the operators Lj accordingly 
to (ff), (nob and operators Bj by formulae 

B-i = 1, Bj = (D x + bi-ij)Bj-i. 

Applying the relations 

TiBj = —ciiBj — HijBj-i mod [L], j — 1, . . . , n 

to the equation Til* = T mod [L] where /* = 5Z?_q Cj-Bj-i one obtains the 
sequence of equations 

(OiTi + lJKn) = 0, 

(oiTj + l)^-) = -^j+i-Hii, i = l,---,n-l, 

^0 = — £l-ffi,0> 

%) = 
which implies contradiction /» = 0. I 

In the totally discrete case the following statement can be proved in similar 
manner. 

Theorem 2 Let equation jj^) possesses the i-integral (j -integral) of the order 
n. Then the Laplace invariant hijk = for some k, < k < n (hij t —k = 0, 
< k < n). 

6 Higher symmetries 

The notion of symmetry for the equations (fil), (0) and (ph is introduced uni- 
formly as any function on the dynamical variables of the given equation which 
satisfies the characteristic equation L{g) = where L is linearization operator. 



11 



Less formally, symmetry can be understand as equation ut — g compatible with 
the given one. It is known Q that equations (||) possessing integrals admit 
also the rich family of the symmetries of the form 

g = P(Y) + Q(X) 

where P and Q some explicitly computable differential operators on D x and D y 
respectively and X, Y are arbitrary x- and y- integrals. 

The aim of this Section is to obtain the analogous formula in the discrete 
case. Let us prove the following Lemma as a preliminary. 

Lemma 3 Let the lattice M) possesses integrals then the kernels of operators 
D x — fui+n Ti — f Ui x are not empty in the space of the functions on dynamical 
variables. Analogously, if the equation M) possesses integrals then the kernels 
of the operators Ti — f Ui +1 , Tj — fi+i.j are not empty. 

Proof. It is sufficient to differentiate the relation defining the integral on the 
higher dynamical variable. For example, differentiating on uf 1 of the relation 
(Ti - l)I{x,Ui, . . . ,u\ n) ) = where u^ l) = D^f) in virtue of the lattice (0) 
yields (f Ui JTi — 1)(7 <„>) = 0, therefore 1/7 <„> belongs to the kernel of the 
operator T, — f Uim . I 

Now let us prove the theorem about the symmetries of the lattice (|l]) (we 
restrict ourselves by this case since the arguments for the equation (0) are quite 
analogous) . 

Theorem 4 Let the lattice (Oj admits an x-integral of the order m and an 
i-integral of the order n. Then the differential operator P and the difference 
operator Q exist of orders not greater than m — 1 and n — 1 respectively such 
that the formula 

g = P(I) + Q(X) 

defines the symmetry of the given lattice for any i-integral L and x-integral X. 

Proof. Accordingly to the Theorem 1, the existence of the x-integral implies 
that hi —j = for some j, < j < m. Applying the Laplace i-transformation 
(|12|), ( fL3| ) for j — 1 times one obtains the general solution of the equation L(g) = 
in the form 

g = -—{D x + b % )-...- — (D x + 6 i+J -_ 2 )(V) 

where ip satisfies equation 

(D x + b i+j - X ){Ti + aij_i)(^) = 0. 



12 



Let (T{ + aij-i)(ip) = 0, then ip can be represented as ip = Kij-o. ■ ■ ■ Kifitpl 
where / is an arbitrary i-integral and <p satisfies equation (T\ + a,i)(ip) = 0. 
Solvability of this equation was proved in Lemma 3 (notice that it is provided 
by i-integral) resulting in particular solution of the form g = P(I) of the char- 
acteristic equation. 

Now consider Laplace ^-transformation. Let the invariant hij vanish, < 
j < n, then the solution of the characteristic equation can be represented as 

g = -!-(2i + on) • . . . • j^—(Ti + Oi){fl>) 

where tp satisfies equation 

{T i + a i )(D x + b i _ lJ )(iP)=0. 

Moreover, accordingly to (|10|), 

h-ij — bi-j-i,o ~ Dx log(J?i_i,_y_i . . • Hi-jfi) 

and therefore the function ip = Hi-\^-\ . . . Hi-jfiLpX is the particular solu- 
tion if X is arbitrary cc-integral and tp satisfies equation ip x + bi—j—ifitp = 0. 
This equation is solvable according to Lemma 3, and one obtains the particular 
solution of the determining equation of the form g = Q(X). I 

As an example we present the formula for the symmetries of the discrete 
Liouville equation (jig): 

u 2 
Ui ,t = D x {uJ) + {Ti - l)((u t + —i-)X). 

Ui+l 

The first term in this formula describes the family of the evolution equations for 
which the lattice ( fig ) defines the Backlund transformation. Note for comparison 
that in the continuous case operators P and Q for the symmetries of the equation 
(logu) X y = 2m are P = D x u, Q — D y u. 

Authors are grateful to A.V. Zhiber, V.V. Sokolov, I.T. Habibullin, S.P. Tsarev 
and A.B. Shabat for interest to the work and useful discussions. The work is 
supported by Russian Foundation for Basic Researches. 
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